Searching for Stable Na-ordered Phases in Single Crystal Samples of 7-Naa:Co02 
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We report on the preparation and characterization of single crystal 7 phase NaxCo02 with 
0.25 < X < 0.84 using a non- aqueous electrochemical chronoamperemetry technique. By care- 
fully mapping the overpotential versus x (for x < 0.84), we find six distinct stable phases with Na 
levels corresponding to x ~ 0.75, 0.71, 0.50, 0.43, 0.33 and 0.25. The composition with x ^ 0.55 
appears to have a critical Na concentration which separates samples with different magnetic behav- 
ior as well as different Na ion diffusion mechanisms. Chemical analysis of an aged crystal reveals 
different Na ion diffusion mechanisms above and below Xc ^ 0.53, where the diffusion process above 
Xc has a diffusion coefficient about five times larger than that below Xc. The series of crystals 
were studied with X-ray diffraction, susceptibility, and transport measurements. The crystal with 
X = 0.5 shows a weak ferromagnetic transition below T = 27 K in addition to the usual transitions 
at T = 51 K and 88 K. The resistivity of the Curie- Weiss metallic Nao.7iCo02 composition has a 
very low residual resistivity, which attests to the high homogeneity of the crystals prepared by this 
improved electrochemical method. Our results on the various stable crystal compositions point to 
the importance of Na ion ordering across the phase diagram. 

PACS numbers: 61.10.Nz, 61.50.Nw, 66.30.-h, 75.20.Hr, 75.30.Cr, 71.27.+a 



I. INTRODUCTION 

The lamellar cobaltates Naa^Co02 have attracted much 
attention since the discovery of superconductivity in 
the hydrated composition Nao.3Co02T.4H20.^ The com- 
plete phase diagram of Naa;Co02 shows a variety of rich 
phenomena, such as A-type antiferromagnetic (AF) or- 
dering for X ~ 0.75, large thermopower enhancement for 
X = 0.85, and AF ordering with a metal-to-insulator 
transition for x = 0.5.^^^^^^^^^ Naa^Co02 is composed of 
alternating layers of Na and C0O2, where the Na ions 
are surrounded by six oxygens which form a prismatic 
cage in the 7 phase. ^ Sodium de-intercalation in pow- 
der samples has been effectively performed via chemi- 
cal extraction, such as a topochemical process using a 
Br2 / acetonitrile solution^ We have demonstrated previ- 
ously an electrochemical method as an alternative route 
to produce homogeneous compounds using a controlled 
overpotential (above the open circuit potential).^ Us- 
ing an aqueous NaOH/H20 solution as the electrolyte, 
high quality single crystals were obtained, and measure- 
ments of their conducting, superconducting and mag- 
netic properties have been reported. ^^^^ However, de- 
tailed studies of the voltage versus Na concentration di- 
agram for X < 1 /2 have been hindered due to difficulties 
related to an oxygen evolution side-reaction in the aque- 
ous solution.^ 



Na ordering - in particular its impact on the geometry 
of the Fermi surface - is a key issue that requires further 
study in order to better understand the transport and 
magnetic properties. ^^^^^^'^^^^^ Local density approxima- 
tion (LDA) calculations indicate that a band with ai^ 
character should create a large hexagonal Fermi surface 
centered around the F point. The existence of such a 
large Fermi surface has been verified from angle-resolved 
photo-emission(ARPES).^^'^^ However, subtle effects due 
to possible superstructure formation and/or strong elec- 
tron correlation remain elusive, such as the missing small 
hole pockets predicted by LDA calculation. Bobroff et al. 
proposed a nesting scenario based on an orthorhombic 
reduced Brillouin zone caused by Na ordering to inter- 
pret the antiferromagnetic and metal-to-insulator transi- 
tions found in Nao.5Co02.^ A reconstructed Fermi sur- 
face may also explain the small hole pocket found in 
Nao.5Co02 and Nao.3Co02 by Shubnikov-de Haas os- 
cillation measurements. ^'^^ Several superstructures have 
been observed and modeled from electron diffraction 
studies of Naa;Co02 with x between 0.15 and 0.75^^ 
Calculations based on electrostatic interaction between 
Na ions alone have been used to propose several types 
of Na ordering patterns. Other calculations point to 
the importance of the coupling between Na- vacancy or- 
dering and Co^+/Co^+ charge ordering, indicating that 
Na+-Na+ repulsion and the Na(l)-Na(2) site energy dif- 
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ference alone is not enough to describe the observed 
super lattices. A recent neutron diffraction study and 
associated numerical simulations have elucidated aspects 
of the Na ordering in single crystals of Naa;Co02 with x 
- 0.75 - 0.92.20 

Much still remains to be understood regarding the 
possible Na ordering patterns across the phase diagram. 
Currently, Na ordering has been studied mostly in sam- 
ples with X = 1/2 and x :^ 3/4. However, reproducibly 
preparing precise stoichiometrics of these phases is not 
straightforward. Many early studies of samples with 
X = 0.7 — 0.9 did not pay much attention to the exact 
ordering of Na ions and often assumed substantial de- 
fect formation due to Na loss at high temperatures Al- 
though Nao.5Co02 shows the highest ordering and phase 
stability among the Naa^Co02 series, phase inhomogene- 
ity and phase separation have been discussed based on 
electron microscopy studies of Nao.5Co02.^^^22 A reliable 
procedure for controlling the Na content and the rate of 
de-intercalation of single crystal samples of Naa;Co02 is 
highly desirable, especially in light of the possibility of 
producing Na ordered phases. 

We find that an electrochemical method which care- 
fully controls the overpotential at the sample surface is 
a particularly reliable route that can approach the real 
equilibrium state and, moreover, is suitable for single 
crystal samples. The topochemical method is an ef- 
fective method to produce powder samples as reported 
previously. 2^ However, the Na level in powders is more 
susceptible to Na loss due to the much shorter diffusion 
path and significantly larger surface-to-bulk ratio. In 
fact, we find that over time, Na diffusion under normal 
storage conditions can affect surface layer stoichiometry 
as deep as 50 jam as we discuss below. Here, we have used 
electrochemistry in a non-aqueous solution to carefully 
prepare a series of homogeneous single crystals. Charac- 
terization of these samples reveals that six distinctively 
stable phases fall in the stoichiometry of x ~ 0.75, 0.71, 
0.50, 0.43, 0.33 and 0.25. Most phases are clearly near 
simple commensurate filling fractions, while 0.71 and 0.43 
require further modeling and simulation. Magnetic sus- 
ceptibility and transport measurements have been per- 
formed on these samples and the results are discussed. 



II. EXPERIMENTAL 

All single crystals in this study were cleaved from 
a crystal grown using the floating-zone method as de- 
scribed previously. ^0^2^ In order to minimize any possi- 
ble Na content differences between the surface and bulk, 
we performed all electrochemical experiments using thin 
crystals with thicknesses as small as 100 /im. The elec- 
trochemical cell was assembled using crystal sample as 
working electrode, Pt as counter electrode, and Ag/AgCl 
with sieved liquid junction as reference electrode. In con- 
trast to our previous work, a non-aqueous electrolyte was 
used: IN NaC104 in propylene carbonate. Single crys- 
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FIG. 1: (Color online) (a) SEM image of an "aged" (- 2.5 
years) crystal of Nao.88Co02 cleaved with the c-axis perpen- 
dicular to the plane The left edge corresponds to the surface 
of the as- grown rod. (b) EPMA analyzed chemical concen- 
tration X as a function of distance from the crystal surface 
towards the interior of the rod. Solid lines are fits with a ID 
diffusion model as described in the text. 

tals of Nao.84Co02 with size near 2x2x0.1 mm^ were 
sandwiched into a spring-loaded cell to make the working 
electrode. A DC voltage Vapp was applied between the 
working and reference electrodes at each applied voltage 
step, spanning the range from -0.25 V to 2.5 V relative 
to the Ag/AgCl reference. Equilibration is achieved by 
waiting until the induced anodic current reaches satura- 
tion and maintained for at least 12 hours. For chemical 
analysis of the stoichiometry, electron probe microanaly- 
sis (EPMA) was performed on freshly cleaved inner sur- 
faces. To avoid electron-beam damage, EPMA data were 
collected using a low beam energy of 15 keV for 1 min du- 
ration on each fresh spot. Structural and bulk magnetic 
data were obtained from X-ray diffraction and SQUID 
magnetometry measurements, respectively. For the X- 
ray diffraction, single crystal samples were oriented, and 
the (OOL) reflections at high angles were used to deter- 
mine the c-axis. In-plane resistivity data of our various 
samples were obtained by performing four terminal DC 
measurements on single crystal samples. The crystals 
were cut to 1mm x 0.5mm x 0.1mm average size. The 
leads were attached with DuPont 4922N silver paint on 
a cleaved ab surface and dried in air, yielding contact 
resistances between 5 — 20Cl. 



III. RESULTS AND ANALYSIS 
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1. Sodium ion self diffusion 

When the as-grown single crystal of Naa^Co02 with 
X > 0.7 is not stored properly within a moisture- free en- 
vironment, a white precipitate builds up on the surface. 
We believe the surface Na ions readily react with the wa- 
ter in the air to form white NaOH precipitate. The va- 
cated sites would induce further Na ion self diffusion. In 
order to better understand the diffusion mechanism of Na 
ions, we have examined the Na distribution profile of an 
"aged" single crystal of Nao.88Co02 that has been stored 
in ambient environment of average humidity of 55% for 
about two and a half years. The original single crys- 
tal was prepared using floating-zone method as reported 
previously The " aged" crystal rod had a layer of white 
precipitate built up on the surface of the boule. We per- 
formed EPMA to examine the Na/Co ion distribution on 
a cleaved surface (a6-plane) which extends radially from 
the core area to the edge of the rod as shown in Fig. 
[TJb). An SEM photograph is shown in Fig. HJa), where 
the crystal surface is defined as the zero position. We find 
that the Na concentration rapidly drops from x o:^ 0.88 
in the interior of the rod to x :^ 0.3 at the surface layer 
within a depth of ~50 /im. The white layer outside the 
crystal consists of a rich concentration of Na (Na/Co > 
1.67) and has been confirmed with X-ray diffraction to be 
NaOH. The Na+ ions, with ionic radius ~ 1.02 A, may 
diffuse through the rectangular pores of threshold radius 
^ 0.91 A formed by the oxygen ions of the NaOe pris- 
matic cage.^^ Since Na+ ions diffuse along the a6-plane 
toward the surface boundary through a vacancy hopping 
mechanism, the observed Na distribution shown in Fig. [1] 
would follow Fick's second law. The diffusion topography 
c(r, t) and the diffusion coefficient Do can be described 

by the diffusion equation ^^^^'^^ = in one di- 

mension. Although it is tempting to fit the whole data 
set with a single diffusion coefficient, we note the c(r, t) 
profile has a pronounced crossover near x ~ 0.5. The 
best fit is obtained under the assumption of two differ- 
ent diffusion coefficients above and below the crossover 
point, as indicated by the solid lines in Fig. [TJ Leaving 
the boundary conditions of c(r, t) at the crossover point 
and the surface edge as free fit parameters, the best fit 
values of Do are ~ 3.49 and 0.76 x 10~^^ cw? / s for dif- 
fusion above and below the crossover point, respectively. 
We note that this is close to the lower limit of Dli for 
Li2;Co02 reported in the literature. The diffusion coef- 
ficients extrapolated from single crystals are less suscep- 
tible to errors due to geometric factors and electrolyte 
permeation, compared to powder samples. 

Interestingly, the observed crossover point is Xc ~ 0.53, 
which is just above half filling. The end point at the sur- 
face boundary is near x ~ 0.30, right below 1/3 filling. 
The fit value of Do for x > 0.53 is nearly five times of 
that from below. This indicates two distinctly different 
types of diffusion species dominated above and below half 
filling, which are most likely to be Na vacancy and Na+ 
ion, respectively. For NaajCo02 with x > 0.5, vacancy or- 
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FIG. 2: (Color online) Na concentration x versus applied 
voltage Vapp and open circuit potential (OCP), where x is 
obtained from EPMA analysis. 

dering is expected as demonstrated by the di- vacancy or 
tri- vacancy ordering patterns proposed by Roger et alS^ 
In the multi-vacancy model, energy is gained through 
Na(2) vacancy cluster formation plus promotion of ad- 
ditional sodium ions from the preferred Na(2) site to 
Na(l) site. The larger diffusion coefficient for vacancy- 
dominated diffusion is also reflected in the lower overpo- 
tential required in electrochemical de-intercalation pro- 
cess, as will be discussed below. The concentration near 
the surface boundary is very close to x = 1/3, which 
is exactly the Na concentration that becomes supercon- 
ducting after fully hydrated.^'' The energy gap between 
X = 1/3 and x = 1/4 appears to be large enough to pre- 
vent additional Na loss from the surface. The crossover 
point near Xc = 0.53 also separates the two main classes 
of magnetic behavior within Naa;Co02, i.e. the Curie- 
Weiss metal above and Pauli paramagnetic metal below4 
The observed different diffusion mechanisms and thus dif- 
ferent Na (vacancy) ordering patterns have shown pro- 
found influence on the magnetic and conducting proper- 
ties of Naa;Co02. 



2. Electrochemical de-intercalation 

The electrochemical de-intercalation method relies 
on the sharp concentration gradient at the sam- 
ple/electrolyte interface. Na ions diffuse to the surface 
layer due to the non-zero overpotential until a new equi- 
librated surface chemical potential is reached. A plot of 
of the Na concentration (determined using EPMA) ver- 
sus the applied electrochemical voltage iVapp) is shown 
in Fig. [2] for our single crystal samples. As the applied 
voltage increases, the Na content decreases similar to pre- 
viously reported results on powder samples. ''^^ The sta- 
bility of the resulting phases is evidenced by the wide 
plateau near regions of dVapp/dx = at x near 0.71, 0.5, 
0.33 and 0.25. 
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FIG. 3: (Color online) Representative Laue photographs for 
(a) the starting crystal with x — 0.84 and (b) a de- intercalated 
sample with the lowest Na content x — 0.25. The six bright- 
est spots have been identified to be {100} with d — 2.83A 
indicated by the small triangle. For x — 0.84, there are 12 
spots surrounding {100} in a ring- like pattern. For x = 0.25, 
orthorhombic superlattice peaks can be identified in addition 
to {100}. 

Crystals prepared with applied voltages between 0.5 
and 0.8 V show nearly constant Na content of x = 0.5. 
In addition, higher voltages near 1.3 - 1.6 V and 2.0 - 
2.5 V produce samples with stable phases of x close to 
0.33 and 0.25 respectively, which had been difficult to 
produce due to the high current side reactions in previ- 
ous studies.^ We also measured the open circuit potential 
(OCP) after each equilibrium phase is reached as shown 
in Fig. [1 The OCP is slightly higher than Yapp for x 
> 0.5, but is significantly lower than Yapp for x < 0.5 
due to the side reaction. Since the OCP primarily re- 
flects the surface energy, the trend of OCP > Vapp for 
X > 0.5 suggests the bulk is equilibrating with a Vapp 
that is lower than the surface energy. The results shown 
in Fig. [2] provide us with new insights into the formation 
of phases below x = 0.5 using Vapp instead of OCP. The 
plateaus in the Vapp — x plot indicate phase coexistence 
regions with dx/dVapp = around ~ 0.5 ,1.3 and 2.0 V. 
The corresponding stables phases are near x = 0.50, 0.33 
and 0.25, which are near the simple fractional fillings 
of 1/2, 1/3, and 1/4. Single crystal quality is demon- 
strated by transmission Laue photography as shown in 
Fig. [3] using synchrotron X-rays (A = 0.4959 A), where 
X = 0.84 is the original as-grown crystal and x = 0.25 
is the most de-intercalated sample. A Na superstruc- 
ture indicated by 12 peaks forming a ring surrounding 
(100) is observed, just like that reported in Ref.^^. An 
orthorhombic superstructure is found in x = 0.25, which 
is different from that observed in x = 0.5 with electron 
diffractionjii. Details of the superstructure analysis will 
be presented separately. 

In addition to the expected phases such as x = 1/2, 
1/3 and 1/4, stable phases of x = 0.43, 0.71 and 0.75 are 
also clearly visible in Fig. [2l Although most early reports 
suggest that x = 2/3 is a preferred equilibrated phase, it 
failed to show up in the Vapp versus x data. This result 
is in agreement with the prediction from Density Func- 
tional Theory (DFT) calculation, which suggests x = 2/3 



is not a ground state due to the lack of a commensurate 
Na ordering pattern. Our data indicate that x = 0.71 
is the most stable phase among 0.75, 0.71 and 0.67, as re- 
flected on the significantly wider applied voltage range of 
the dx/dVapp = plateau. The stable phase of x = 0.71 
can be prepared with Vapp between -0.1 - O.IV, while 
X = 0.75 is barely defined by a plateau of width less 
than 60 mV. The x = 0.75 phase has been shown to be 
at the lower boundary of the range (0.75 < x < 0.85) 
for compositions that possess A-type AF ordering at low 
temperature. ^^'^^ In addition, x = 0.75 falls at the ex- 
act boundary that minor x fluctuations would separate 
Na(2) to occupy either high symmetry 2c (2/3, 1/3, 1/4) 
or low symmetry 6h (2x, x, 1/4) sites within the space 
group P63/mmc.^^ Such behavior indicates that x = 0.75 
is likely to be a metastable transient phase rather than 
a truly equilibrated ground state. The x = 0.71 phase 
is very close to a simple fractional filling of 5/7 that has 
been predicted to be a ground state from Density Func- 
tion Theory (DFT) model calculation;^^ however, there 
are many other choices for a simple fraction that are near 
0.71. Roger et al. have proposed an ordering of multi- 
vacancy clusters within the Na layer for x > 0.7, i.e., Na 
ions shift from the preferred Na(2) site to the unfavor- 
able Na(l) site that is directly on top of the Co ions to 
further reduce stabilization energy.^^ 

Another equilibrated phase below half filling is x = 
0.43. Although the dx/dVapp = plateau for x = 0.43 
is less pronounced than that for x = 0.75, the narrow 
plateau sits between two distinctly different slopes as 
shown in Fig. [2l Considering the dominant diffusion 
species for x > 0.5 are vacancy clusters, the diffusion 
species for x = 0.43 would be Na ions. The fact that 
0.43 being close to a simple fractional filling of 3/7 sug- 
gests that Na ions partially fill a superlattice consisting 
of 7 Na ions per unit supercell. Indeed Na trimer order- 
ing of simple hexagonal unit a' = >/7a has recently been 
observed on the surface of an in-situ cleaved Nao.84Co02 
crystal using scanning tunneling microscopy, where half 
of the Na ions (x ~ ^^) are supposed to occupy each 
cleaved surface pair equally.^^ Whether such 3/7 super- 
lattice ordering exists within the bulk requires further 
Laue transmission diffraction analysis. 



3. Phase homogeneity 

One of the major concerns in the experimental study 
of Naa:Co02 compounds is the question of phase ho- 
mogeneity. When Na ions are de-intercalated through 
topochemical or electrochemical routes, the remaining Na 
ions should rearrange to find the ground state for each 
Na layer while maintaining the original CoOe stacking. 
The c-lattice parameters for our series of crystals were 
determined with X-ray diffraction and are summarized 
in Fig. m Since Naa^Co02 samples have highly oriented 
at (OOL) planes even after the crystal is pulverized, X- 
ray diffraction has been performed on aligned crystals 
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FIG. 4: (Color online) The c- lattice parameter of Naa;Co02 
plotted against x. Crystal homogeneity is confirmed by the 
FWHM in 20 of the (008) Bragg peak. The instrumental 
resolution is 0.05 degrees. 



measuring the (OOL) Bragg peaks. The c-axis is known 
to strongly depend on Na content as reported by several 
other groups. ^'^^^^^ Our c-axis lattice size is extracted 
from the (008) peak angle, and our results are consis- 
tent with those reported previously.^ We notice the c- 
lattice parameter depends almost linearly on x with a 
weak deflection near x ~ 0.55. This behavior may sug- 
gest two regimes which follow the empirical Vegard's rule 
for substitutional impurities. The c-axis expansion as a 
function of lower x has also been reflected on the soft- 
ening Ai^ mode along the c-direction^^ As shown in the 
inset of Fig. [H phase homogeneity is evidenced by a con- 
sistently narrow full width at half-maximum (FWHM) 
for 20 of the (008) peak. This width is comparable to 
our instrumental resolution limit, which is determined 
using a Si standard to be 0.05 degrees. The possibility of 
mixed-phase crystals is ruled out within our X-ray reso- 
lution limit for nearly all of the samples prepared by our 
method, except for samples prepared using Vapp between 
0.82 - 0.88 V as shown in Fig. [5l Two-phase behavior is 
observed within the narrow voltage range between 0.82 
- 0.86 V, which corresponds to Na levels slightly below 
X = 0.5. The sample evolves from single phase of x = 0.5 
to a single phase of x = 0.43 through mixed-phase inter- 
mediate states. In contrast, there is no two-phase signa- 
ture observed between single phase x = 0.55 and x = 0.5. 
No other samples exhibited two-phase signatures regard- 
less of whether Vapp falls in the dx/dVapp = plateau 
range. This suggests that phases of solid-solution with- 
out specific Na ordering do exist, at least in the time 
frame of days during preparation and characterization. 
Since the FWHM of (OOL) peaks provide us only with 
coherence information along the c-axis, it is possible that 
in-plane ordering may yield potential microscopic phase 
separated domains as commonly found in this series i^^^ 
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FIG. 5: (Color online) X-ray diffraction scans of the (008) 
peak for different applied voltages near x = 0.5. For x > 0.5, 
all of our crystals appear to be single phase. For x < 0.5, two- 
phase coexistence is observed for the range between x = 0.5 
and X — 0.43 



4- Magnetic properties 

Additional evidence of phase purity is provided by 
magnetic susceptibility measurements. Figure [6] displays 
the complete susceptibility data measured under a mag- 
netic field of 1 Tesla applied along the a6-direction. These 
results are consistent with most of the studies reported so 
far: magnetic behavior evolves from Curie- Weiss-like to 
Pauli-like as Na ions are de-intercalated from 0.84 to 0.25, 
without any trace of CoO or C03O4 impurity^ i^^i^^i^^ 
Assuming identical Van Vleck susceptibility contribution 
for the Pauli paramagnetic metal region, the tempera- 
ture independent portion of the average susceptibility 
(x(300K)) should be proportional directly to the den- 
sity of states at the Fermi level. The susceptibility at 
T = 300 K x(300K) for x < 0.6 are shown in the inset 
of Fig. [6l where x(300K) for x = 0.33 appears as a local 
maximum and x = 0.5 sits at the local minimum. This 
trend agrees well with 7 values extracted from specific 
heat measurement, where 7 is directly proportional to 
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FIG. 6: (Color online) (a) Magnetic susceptibility versus tem- 
perature for crystals measured under applied field of 1 Tesla 
along the a6-direction. The value of the susceptibility at 
T = 300 K has been plotted in the inset, which has a maxi- 
mum near x — 0.33 and a minimum near x — 0.50. (b) The 
slope of the average susceptibility versus T at 150 K is plotted 
against x. There appears to be a crossover near x = 0.55. 

the density of states at Fermi level. The enhancement 
of density of state at Fermi level for x = 0.33 reflects 
its unique superconducting properties after hydration.^ 
The local minimum of x(300K) for x = 0.5 also indi- 
cates its proximity to the insulating state that sets in 
below T = 51 K.^ The slope of the spin susceptibility 
at 150K in Fig. [6^ is an indicator that separates Curie- 
Weiss (negative) and Pauli (positive) behaviors. The 
crossover point is close to x = 0.55 as shown in Fig. [6)3, 
which is consistent with the crossover points that have 
been observed in Fig. [1] of Na self diffusion and c-lattice 
parameter as shown in Fig. [H In Fig. [6^, the in-place 
susceptibility Xa6(T) for x = 0.71 and x = 0.84 are nearly 
identical above ~ lOOK, but x = 0.71 does not show any 
magnetic ordering down to 1.8 K. This behavior suggests 
X = 0.71 is the first equilibrated phase below x = 0.85 
that does not possess the well known A-type AF order- 
ing below 22-28 Magnetic properties of the newly 
identified Nao.7iCo02 and Nao.43Co02 compositions are 
summarized in Fig. The Curie- Weiss-like low temper- 
ature upturn for X = 0.71 agrees very well with the recent 
2^Na NMR Knight shift data on samples with x ^ 0.7.^^ 
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FIG. 7: (Color online) Magnetic susceptibility versus temper- 
ature for Nao.7iCo02 and Nao.43Co02, with H — 1 Tesla for 
both H/ /c and H/ 1 ah directions. The inset shows a low field 
measurement for Nao.43Co02 using 100 Oe under zero- field- 
cooled and field-cooled conditions. 

Powder average of Xah and Xc can be fitted to a simple 
relationship of x(^) = Xo + C/(T — 0)*^ Fitting over 
the temperature range T = 50 - 250 K yields C = 0.078 
cm^K/mole and O ~ -42K. In a simple ionic picture, this 
corresponds to about a 72% concentration of Co^+ ions 
that participate in the weak antiferromagnetically corre- 
lated paramagnetic behavior. Then the electrons on the 
remaining ~28% of the Co ions are itinerant. We specu- 
late that a special Na ion ordering of x = 0.71 may form 
a large supercell that provides the underlying potential 
for such partial charge delocalization. 

Nao.43Co02 shows Pauli-like behavior just like all 
other samples with x < 0.55. However, low field measure- 
ment revealed a weak thermal hysteretic behavior below 
~40 K as shown in the inset of Fig. [71 The hysteresis 
onset near 40K for x = 0.43 is very close to the metal- 
to-insulator transition found in x = 0.5, which seems to 
suggest the inclusion of nearby x = 0.5 phase. In addi- 
tion, there are traces of weak 51 K and 88 K anomalies 
observable at the level of a few percent. However, the 
transition at ~51 K for x = 0.5 does not show thermal 
hysteresis, instead, the onset of hysteresis is found to set 
in below the weak ferromagnetic phase ~27KJ^ Based 
on the room temperature structural evidence shown in 
Fig. [5l we can produce samples with no x = 0.5 phase 
inclusion within the x = 0.43 matrix. We cannot rule 
out the possibility of phase segregation below room tem- 
perature. The observed thermal hysteresis of the sus- 
ceptibility below ~40 K may be related to either weak 
ferromagnetism or spin glass behavior. For Nao.5Co02, 
three phase transitions can be identified at ~27, 51 and 
88 K in the magnetic susceptibility measured along the 
a6-direction under low magnetic field, as shown in Fig. [HI 
The transition at 88 K is due to a novel antiferrmagnetic 
ordering transition, where the ordered spin arrangement 
and direction has been deduced from polarized neutron 
diffraction.^ The transition at 51 K correlates to a metal- 
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FIG. 8: (Color online) Magnetic susceptibility versus temper- 
ature for Nao.55Co02 and Nao.5Co02 prepared using 0.44 V 
and 0.50 V for 24 hours each. The applied field is 100 Oe 
along the a6-direction, and measurements were taken under 
zero- field- cooled and field-cooled conditions. The inset shows 
similar measurements with a higher field of 1 Tesla. 

to-insulator transition which may be related to nesting 
of the Fermi surface.^ Most of the magnetic measure- 
ment results of Nao.5Co02 at high field before show low 
temperature upturns below ~30 K of various sizes, usu- 
ally attributed to trace amount of isolated paramagnetic 
spins 4^1^ Interestingly, we find the transition at 27 K 
shows ferromagnetic character with thermal hysteresis at 
low field as well as field hysteresis at 5 K (not shown). 
The significant thermal hysteresis found below ~27K for 
Nao.5Co02 strongly implies its weak ferromagnetic na- 
ture. This transition below 27 K appears to further in- 
crease the in-plane resistivity^^ 

The Nao.55Co02 composition has been previously re- 
ported to have a strong in-plane ferromagnetic ordering 
below ~20K with no magnetic anomalies found near 51 
or 88 K.^^ One immediate concern on preparing pure 
NaajCo02 crystal with x = 1/2 is how to separate 
the potentially mixed phases between Nao.55Co02 and 
Nao.5Co02. A series of fully equilibrated crystals have 
been prepared with different applied voltage in small 
steps between 0.40- 0.50 V; the results of structure evo- 
lution are shown in Fig. [5] and magnetic susceptibilities 
are plotted together with that of x = 0.5 in Fig. [HI High 
field susceptibility measurements show a clear difference: 
nearly temperature independent behavior for x = 0.55 to 
the one with positive slope for x = 0.5. Low field mea- 
surements indicates our Nao.55Co02 crystal has a weak 
ferromagnetic transition near 27 K without any 51 K or 
88 K anomalies, in contrast to Nao.5Co02. This rules 
out the possibility of a mixed phase within such a small 
concentration range. Nao.55Co02 and Nao.5Co02 should 
be considered as two distinctly different phases, although 
the former is closer to a critical transient state with an 
extremely narrow dx/dVapp plateau as indicated in Fig. 
[2j We note that the saturated moment and coercive 
field for Nao.55Co02 obtained from the electrochemical 
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FIG. 9: (Color online) (a) Temperature dependence of the 
in-plane resistivity for several of the NaxCo02 crystals in the 
series, (b) Low temperature resistivity of the composition 
with X = 0.71. Two regimes are clearly seen: a tempera- 
ture dependence below 5 K (inset) and a T-linear temperature 
dependence between ~ 5 K and 80 K. 



method is lower than that prepared by the chemical ex- 
traction method using iodine/acetonitrile solution. In 
Nao.5Co02, the size of the magnetic susceptibility upturn 
below ~27 K appears to be highly sample dependent 
in the literature. ^^^^^^^^^ One major difference between 
chemical and electrochemical methods has been com- 
pared on LiajCo02. Chemical de-intercalation method 
has a significantly higher charging rate, and thus easier 
to access the met ast able states. We believe the chemi- 
cal method with a much higher charging rate may yield 
metastable phases of Nao.55Co02 of higher weak FM mo- 
ment, and these may eventually fall into a more equili- 
brated phase similar to that obtained using the electro- 
chemical route. 
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5. Transport properties 

The transport properties of four of our newly synthe- 
sized compositions are summarized in Fig. [9l As shown 
in Fig. [9](a), the compositions with x = 0.71,0.55, and 
0.25 exhibit metahic conductivity, and are consistent 
with the trend suggested by the phase diagram previ- 
ously proposed.^ The sample with x = 0.43, does not 
behave like a good metal and appears to have a large 
residual resistivity. This particular sample may contain 
coexisting regions of different phases, such as x = 0.5, 
which are more stable at the local scale. Indeed, after 
annealing at = 400K for several hours (not shown), 
this sample displays a transition to an insulating regime 
below T ^ 50K, which is a particular feature of the 
X = 0.5 composition^. The resistivity of the x = 0.71 
sample at lower temperatures is shown in detail in 
Fig.[9]^b). A roughly T-linear temperature dependence is 
observed between ~ 5 K and 80 K, consistent with pre- 
vious measurements.^ Fermi liquid behavior is recovered 
at low temperatures where electron-electron scattering 
dominates, as noted previously. Here, below ~ 4K the 
data may be fit by the function of p{T) = po + AT^. 
Our fit (plotted as the solid line) yields values of 
p = 15.61 pn - cm and A = 0.2226 pQ - cm K''^. 
The residual resistivity po may arise from scattering 
due to domain boundaries and/or perturbations to the 
periodic potential felt by electrons within the ah plane, 
such as defects in the Na arrangement. The observed 
value of po is about three times smaller than previous 
measurement s"^"^ and attests to the high homogeneity of 
the current Na-ordered x = 0.71 sample. The value of 
A, while somewhat smaller than previous measurements, 
is still relatively large and indicates an large Kadowaki- 
Woods ratio, kxw = Aj^^^ assuming 7 is similar to 
previous measurements^ 



IV. SUMMARY 



In summary, we have demonstrated that a non-aqueous 
electrochemical de-intercalation process can be used to 
prepare high quality single crystal samples of 7-Naa;Co02 
with well-defined Na superlattices. The sodium ion self 
diffusion mechanism has been demonstrated to be well fit 
by a ID diffusion model on an aged single crystal. The 
dominant diffusion species above Xc ^ 0.53 is consistent 
with Na vacancies, and the diffusion coefficient is higher 
for X > Xc, compared to x < Xc- The existence of sta- 
ble phases are reflected in the dx/dVapp = plateaus, 
among which x = 1/2, 1/3 and 1/4 are the most sta- 
ble. The data suggest that x ~ 0.75, 0.55 and 0.43 are 
somewhat less stable. On the other hand, x ~ 0.71 is a 
very stable phase with a sizeable dx/dVapp = plateau 
width. This stable x = 0.71 phase is intriguing for its 
potential Na ordering but without any magnetic order- 
ing down to 1.8K. The absence of A- type AF ordering 
for X = 0.71 and the presence of weak FM in x = 0.5 
and 0.55 may be related to subtle interactions between 
the Co ions and the specific potential caused by Na ion 
ordering, as it is known that Na superstructures can af- 
fect the Fermi surface geometry. Clearly more work is 
required to investigate the effects of Na order across the 
Naa^Co02 phase diagram. Our progress with this electro- 
chemical technique opens the door for further studies of 
the physical properties of the stable Na ordered phases 
in single crystal samples, 
"f" corresponding author: fcchou@ntu.edu. tw 
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